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DEVELOPMENT OF A SIMPLIFIED PROCEDURE 
FOR THRUST CHAMBER LIFE PREDICTION 

SUMMARY 


An analytical design procedure for predicting thrust 
ch€unber life considering cyclically induced thinning and bulg- 
ing cf the hot-gas-wall is developed. .The hot-gas-wall, com- 
posed of ligaments connecting adjacent cooling channel ribs and 
separating the coolant flow from the combustion gas, is sub- 
jected to pressure loading and severe thermal cycling. Thermal 
transients during start-up and shut-down cause plastic straining 
through the ligaments. The primary bending stress superimposed 
on the alternate in-plane cyclic straining causes incremental 
bulging of the ligaments during each firing cycle. This basic 
mechanism of plastic ratcheting is analyzed and a method 
developed for determining ligament deformation and strain. 

The method uses a yield surface for combined bending and mem- 
brane loading to determine the incremental permanent deflection 
and progressive thinning near the center of the ligaments which 
cause the geometry of the ligaments to change as the incremental 
strains accumulate. Fatigue and tensile instability are af- 
fected by these local geometry changes. Both are analyzed and 
a failure criterion developed. Results of the simplified anal- 
yses are shown to compare favorably with experimental data and 
finite element analysis results for OFHC (Oxygen Free High 
Conductivity) copper. They are also in reasonably good agree- 
ment with experimental data for NARloy Z, a copper- zirconium- 
silver alloy developed by the Rocketdyne Division of Rockwell 
International . 



DEVELOPMENT OF A SIMPLIFIED PROCEDURE 
FOR THRUST CHAMBER LIFE i'REDICTION 

1.0 INTRODUCTION 

Life predictions of regeneratively liquid cooled rocket 
thrust chambers have been based on low cycle fatigue principles. 
Tests of thrust chambers, however, [1, 2]* revealed that coolant 
channel walls in the failed areas have exhibited progressive 
incremental thinning and bulging during the heating and cooling 
cycle associated with each firing. Failure analyses performed 
as a part of this work indicate that ductile rupture is a more 
limiting mode of failure than fatigue failure. Thus, material 
tensile instability was considered in addition to fatigue 
damage mechanisms in the present structural evaluation. 

The mechanism of hot-gas-wall or ligament incremental 
distortion can be investigated by inelastic finite element analy- 
ses. For the cyclic histories of interest, however, such analy- 
ses are difficult to perform and require extensive computer 
time. Their use for a single characteristic case is feasible 
and provides valuable information. In fact, such analyses 
were carried out within NASA-Lewis research programs for simpli- 
fied histories of thermal loading. Such methods, however, are 
not suitable for evaluating the effects of changing various 
design and operating parameters such as geometric configura- 
tions, material properties, pressures and temperatures. More 
general methods which do not require lengthy inelastic finite 
element solutions are therefore needed. 

The present work is aimed at; (1) analyzing the failure 
mode of the cyclically loaded thrust chamber, (2) developing 
simplified but conservative methods for evaluating the strains 

♦Numbers in brackets designate references at end of report. 
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and deformations, and (3) providing methods for predicting the 
cyclic life of the thrust chamber. 

One or two cycles of elastic and inelastic analysis were 
analyzed to determine stress and strain resultants for the 
simplified models. Other available inelastic finite element 
results and test data were also used to determine the accuracy 
of the simplified models. 

The scope of the work is as follows: 

1. Develop a simplified analytical model for 
determining hot-gas-wall deformation and strain. 

2. Determine accuracy of simplified model by 
comparison with inelastic finite element 
analysis results. 

3. Develop a failure mode evaluation criterion. 

4. Provide a simple design procedure for making 
life predictions for thrust chambers. 

Each of these items will be discussed in the sections to ^‘>llow. 
Also the assumptions made in the development of the analysis 
procedure for each item are disr iissed in their respective 
sections . 

2.0 PLUG NOZZLE THRUST CHAMBER 

The plug nozzle thrust c.iamber shown in Figure 1 is con- 
sidered. Figure 1 provides dimensions for the thrust chamber 
tested at NASA Lewis Research Center. The sample solution and 
computer finite element calculations given in the J^pendices 
were performed for this thrust chamber. However, the resulting 
simplified method of analysis is applicable for a general case 
where the geometry, loading conditions and materials can be 
altered to comply with the design of a specific engine. The 
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Thrust Chamber Assatnbly 



152.4 nun 
(6.00 in.) 



rvlinder Cross Section 


Analytical Model 


4.953 mm 
(.195 in.) 


1.27 mm 
(.05 in.) 




Radially oriented 
Mid-Plane of Coolant Channel 
Closeout Wall 


i6=1.6876nin 
(.066 in.) (1-30 in.) 

2It®.889 mm 
(.035 in.) / 

Combustion Gas £_ inner Wall with n»72 
Coolant Channels 


FIGURE 1 PLUG NOZZLE THRUST CHAMBER 
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plug nozzle assembly consisting of the contoured centerbody and 
flanged cylinder is shown along with cross sectional details of 
the cylinder. The inner liner of the cylinder contains ?xial 
flow coolant channels of constant cross section and is con- 
structed from either OFHC copper, half -hard Amzirc (American 
Metal Climax, Inc. copper alloy) or NARloy-Z. The closeout wall 
is electroformed copper in all cases. Since no failures occurred 
in the centerbody only the cylindrical part of the nozzle was 
investigated. 

Typical firing cycle history data for the thrust ch£unber 
is illustrated in Figure 2. The data refer to the critical 
section where failures were observed during testing. Only this 
section was considered in the analysis. The first cycle begins 
with a cooling period when the liquid hydrogen enters the chan- 
nels. The uniform temperature 29**K (53"R) is reached before a 
sudden rise of temperature beginning at the ignition point of 
1.5 seconds. 

The mechanical properties used to characterize the cylinder 
wall materials were provided by NASA-Levis. Reference 3 provides 
temperature dependent thermal expan&xon, modulus of elasticity 
and static stress-strain properties of half-hard Amzirc, NARLoy-Z, 
annealed OFHC and EFCU (electroformed copper) . 

3.0 DEVELOPMENT OF A SIMPLIFIED ANALYTICAL MODEL FOR CALCULATING 
HOT-GAS-WALL DEFORMATION AND STRAIN 

3.1 Simplified Model 

The inner wall shown in Figure 1 separating the coolant 
flow and the combustion gas is subjected x .3 severe thermal 
cycling. The temperature difference between the closeout 
and inner walls that occurt during the thermal transient 
causes plastic straining through the ligament within the 
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FIGURE 2 THE TIME HISTORIES OF HEAT TRANSFER COEFFICIENTS 
AND TEMPEivU.TURES FOR A CHILLOOWN AND FIRING CYCLE 



range of 1 to 2 percent. The temperature history for the 
closeout and inner walls at the middle plane of the cooling 
channel is shown for one cycle in Figure 3. The inner wall 
is attached to the closeout wall by radial ribs. The outer 
wall remains elastic within the entire thermal cycle forc- 
ing the inner wall ligaments into captive plastic straining. 
Straining in the hoop direction is accompanied by high 
axial strains. 

The ligaments separating the coolant and the combus- 
tion gas are subjected to coolant pressure load and ther- 
mally induced through-the-wall bending. The resulting 
stresses interact during plastic straining causing incremen- 
tal bulging or ratcheting directed radially inward. The 
mechanism of pressure induced ratcheting is illustrated in 
Figure 4 . Primary bending stresses persist during alter- 
nate in-plane cyclic straining resulting in a small repeti- 
tive incremental permanent deflection of the ligaments 
during each firing cycle. The simple model in Figure 4 
does not include interaction of shear stress due to pressure 
and the thermally induced bending. T^?ir effects are evalu- 
ated later in the report. 

3.2 Analysis of Basic Mechemism of Plastic Ratcheting 

The approach used to determine increments of plastic 
strain within the load cycles is explained in Figure 5. 

The axial forces and bending moments acting on the rectan- 
gular cross section as shown in Figure Sa are considered. 

The bar simulates the response of the ligament of a cooling 
channel between the two ribs. The bar is cyclically 
stretched and squeezed plastically in the presence of sus- 
tained bending. Note that the "X" direction in the beam 
model corresponds to the hoop direction in the ligament. 
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(Axial Strain) 


FIGURE 4 ELASTIC RATCHETING DUE TO INTERACTION 
OF PRESSURE AND THERMAL CYCLING 
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X - Hoop Direction 
y - Axial Direction 



FIGURE 5 FUSTIC RATCHET INCREMENTS IN A BAR 0/ RT'CTAN 
SULAR CROSS SECTION SUBJECTED TO SUSTAINED 
BENDING MOMENT AND CYCLIC AXIAL STRAIN 

a) Ligaioent Bean Element Loading 

b) Interaction Curve 
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To obtain a conservative bound on the plastic strain 
incresents, strain hardening of the Material is ignored. 

The yield curve for a beam subjected to bending and 
axial force is composed of two parabolas as shown in Fig- 
ure 5b emd is defined by the relation* [4] : 

F = iB + n^-l = 0 (1) 

idiere n and n are dimensionless vauriables defined by: 

m = i!- (2a) 

o 

n * (2b) 

"o 

In the abovs* M and N denote the bending moment and axial 

force* %diile and N denote the yield bending moment 
o o 

and yieZd axial force given by: 


N 


O 



(3a) 


N « 2HS 
o y 


(3b) 


for a rectangular be£un of unit width* height 2R and yield 
stress Sy. 

Plastic flow vectors q are normal to the yield curve. 
The curve is plotted using generalized stresses n and n 
related to axial 'and bending tractions* respectively. The 
associated generalized strain rates A and 0 are defined by 
the average hoop strain rate e and the rate of curvature of 
the bent bar* K* respectively* as shown in Figure 5. If 
hoop increments of plastic strain are Icnown the correspc 
ing increment of curvature is defined by the slope of the 
yield surface and can be determined (based on the normality 
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of resulting strain rate vector g) as shown in Figure 5b. 

The generalized strain rates can be obtained by 
partial differentiation of equation (1) as follows: 

i = a|^ = 2nA (4a) 

e = a|| = A (4b) 

where A is an arbitrary positive scalar. Thus, for the 
interaction of hoop and bending loads, the relation between 
X and 6 components of the q vector is simply: 



For known values of bending m %diich remain constant 
within the cycle, n can be obtained from the yield line of 
equation (1) . Hoop increments of plastic strain are being 
reversed within t«fo halves of the thermal cycle. However, 
plastic increments of curvature within each half of the 
thermal cycle are of the same sign. Thus, they accumulate 
causing incremental bulging of the cooling channel wall. 

3.3 Interaction of Shear Stress 

For high pressures in the cooling channels the model 
of the .thrust chamber must also include shear stress t . 

The method of solution is analogous. However, the solu- 
tion includes generalized shear stress and the correspond- 
ing sheSr strain. 

The yield surface based on Tresca's yield criterion 
is given by Peterson tt at. 

F«m + 2^-k»0 (6a) 
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where 


and 



(6b) 


s - ^ is the dimensionless shear stress. 


The generalized strains X , 9 and # are related to 
the hoop strain, curvature and shear strain, respectively. 
Using the normality law: 


\ 


e 




♦ 




(7a) 

(7b) 

(7c) 


The proportions between the strain con^nents are thus; 

3F , IF ^ IF , : (2n) : f 2s - (8) 

am an as \y / ^ ^ _ g2 j 

The generalized variables X, 9 and 4 can be related 
to the hoop strain, curvature and shear strain by consid- 
ering the rate of dissipation D given by [4] ; 

D * M K + N e + S Y (9) 

where M, N and S denote the bending moment, hoop force 

and shear force, respectively, and K, e and y denote the 
curvature rate, hoop strain rate and shear strain rate, 
respectively. 
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Using equation (2) , equation (9) can be written as: 


0 s mM K + nN c + sS Y 
o o o 


where 


s = 




(10a) 


(10b) 


since the cross sectional area for vinit width equals 2H. 


Then using (3), (10a) and (10b): 


D = Sy[m(H2K) +n(2He) + s(Hy) 1 (11a) 

Also in terms of the generalized strain rates, the rate 
of dissipation is given by: 

D = Sy(m e + n X + s ♦] (lib) 

Thus the relationship between the generalized strain rates 
and the hoop strain rate, curvature rate and shecu: strain 
rate are, respectively: 


X = 2HC 

(12a) 

0 = H^K 

(12b) 

♦ = hy 

(12c) 


For known increments of hoop strain X and known 
values of m and s, the remaining components of strain 6 
and I) (shear) are determined by the ratios (8) . Since the 
load point must remain on the yield surface, the known 
values of m and s determine n which is obtained from 
equation (6a) . 
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The curvature and shear strain as determined from (12) 
can then be integrated along the length of the ligament to 
obtain the bending and shear deflections. Once the total 
deflection is known, the corresponding thinning in the 
ligament is calculated by assuming that the volume of the 
material remains constant. Integration of the deflection 
and determination of ligament thinning is described in 
Section 3.4. 

It can be shown from the generalized yield surfaces 
derived in Ref. [5] that though the yield surface given by 
equation (6a) is derived for a uniaxial model, it also holds 
for biaxial conditions when ri 2 < n^^ (n^^ ® n * hoop stress 
variable in equation (6a)) and m and s are small, as in the 
present case. Thus the generalized strain relationships 
derived ad>ove are also applicable for the biaxial case. 


3.4 Deformation History 

From equations (7a) , (7b) and (12a) , (12b) : 

1 » 1 * V 1 - H K ^ H K 


For the complete cycle, the plastic strain generated by 
the tenqperature difference during heat-up in the hoop 
direction is fully reversed when the temperature becomes 
uniform at the end of the cycle. Thus.- 



v’here ACp is given by sum of plastic strain range due to 
differential thermal expansion Ael and the correction due 


to thermally induced bending Aep^ which are evaluated in 
Sections i.5 and 3.6, respectively. 


Thus, the curvature is given by: 
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( 15 ) 



Sindliarly from equations (7a) , (7c) and (12a) , (12c) : 

f = /2s - -T-^- - ] /2n = V2e = Y/2e (16) 

X X ^ v/l-s2 j 


Making use of equation (14), the shear strain is then: 



(17) 


The curvature and shear strain as determined from 
eqiiations (15) and (17) for each cross section can then be 
integrated along the length of the ligament to obtain the 
bending and shear deflections. Once the total deflection 
is known, the corresponding thinning in the ligaunent is 
calculated by assinning that the volume of the material 
remains constant. 

Experimental evidence shows that the deformed shape 
of the ligament can be approximated by a linear variation 
in thickness as shown in Figure 6. If 5 denotes the deflec- 
tion per cycle, the thinning after N cycles is then given 
by: 


. ^ N g w 

^ (i + w) 


(18) 
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Note that even though the incremental deflection remains 
constant, the incremental strain at the center of the liga- 
ment increases with each cycle due to ligament thinning. 




3.5 Stress and Strain Resultants 


The analysis of incremental plastic deformation in- 
volves generalized stresses m, n and s and the associated 
generalized strains. The generalized moment and shear and 
the generalized strain in the hoop direction are needed as 
input data tor the solution. Their numerical values can 
be obtained by crude hand calculauions using results of 
the thermal analysis and elastic finite element calcula- 
tion or by performing the full inelastic analysis of the 
first two load pycles . 


Hand Calculations 

The bending moment distribution along the beam can be 
approximated by the solution for a clamped beam under 
uniformly distributed load: 


M 




(19) 


where p is the unit load, I is the length of the beeun and 
X is the axial coordinate measured from the rib as shown in 
Figure A-1 of Appendix A. 

Shear force is simply: 


S 



( 20 ) 


The average hoop strain range, Ae^ , in the ligament 
i.'- evaluated based on the differential thermal expansion 
of the ligament at average temperature Tj^ and the closeout 
wall at average temperature T^. Since the thickness of the 
'igament is small with respect to the thickness of the 
closeout wall, the elastic deformation of the closeout wall 
can be disregarded and the range of plastic strain in the 
hoop direction is: 
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Ae’ 

Pi 


(T. a. - T ot ) - (T. a . - T a ) . 

11 o o max 11 o o min 


.]■(' 


+ s 


max 


min 


j 


/E 


(21) 


and are ♦;he thermal expansion coefficients of the 


where 

ligeunent 6uid closeout., respectively, and 


(T.o. - To) and 
11 o o max 

^^i°*i ” ^o°‘o^min maximum and minimum thermal strains, 


max 


respectively, that occur during the loading cycle. 
and are respectively the ligament material absolute 

yield strengths corresponding to the ligament average tempera- 
tures at the times in the cycle when 




o o max 


(T.o. - are calculated. 

11 o o min 

The average temperature uf the ligament can be assess- 
ed assuming steady heat transfer conditions between the 
hot gas at adiabatic temperature and the evaporating hydro- 
gen in the cooling channel. 


Using a one dimensional heat flow model ad assuming 
that heat enters the metal through the inner surface of 
the thrust chamber and is being received by two sides of 
the cooling channel as shown in Figure 7 , the resulting 
temperature of the sides of the ligament are: 


/4H 


T. = 


A, + A, 


'*'^2 


) W 


1 + T4 


1 + 


••lAi 




4H 1 

k Ajj + Aj + A 2 


h,A, + h A 

11 2 2 




(22a) 


T, * T, + 


h.A +h.A 


.( 


4H 1 1 

k Ajj+A^ +Aj hjAj+h^A^ 




+ T. 


1 1 


2 2 


1 + 


ha r^H 


X 


+ A, +A, 


1 


(22b) 

where h and A denote heat transfer coefficients and sur- 
face areas as shown in Figure 7. 
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One Dimensional Heat Flow Model ; 

. ^ A1+A2 

hfl Ajj(Ti-T2) -2—^ (T2-Ts) = (hiAi+h2A2)(Ts-T.,) 


FIGURE 7 THERMAL MODEL OF LIGAMENT 
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Thus the average temperature of the ligament is: 


Ti = 


T2 + T3 


( 23 ) 


A certain part of heat bypasses the cooling channel and 
heats the closeout wall, thus reducing the differential 
straining. The temperature of the closeout wall is approxi- 
mately proportional to the average temperature of the 
ligeunent: 

T = AT. (24) 

o 1 


The approximate value of constant A for the considered 
thrust chamber is 0.35. This value for A was obtained 
for a specific geometry and will probably be different 
for other geometries. 

The temperatures T 2 and T^ can also be used to deter 
mine the thermally induced bending stress through the 
ligaments which for equibiaxial stress field is: 


ATaE 

^b ~ 2(1 - v) 


(25) 


where AT is the temperature difference through the ligeunent. 


Finite Element Thermal and Elastic Analyses 

A more accurate evaluation of stress and strain 
resultants may be obtained by performing thermal and stress 
emalyses for the extremal conditions using the finite ele- 
ment method. The formulation of the problem and instruc- 
tions to prepare the input data for the MARC program are 
given in Appendix B. Selection of thermal steps is 
optimized by the program. The stress run should be per- 
formed for the instants resulting in the maximum range of 
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stress. The calculation of bending and shear tractions is 
straightforward. The plastic strain range in the hoop 
direction should be obtained by averaging the total 
strains through the thickness and correcting the result 
to account for elastic unloading as follows: 




( 26 ) 


In order to obtain conservative evaluation of plastic 
strain the correction for unloading derived by the second 
term in equation ( 26 ) is given for the equibiaxial state 
of stress. 

The thermal bending through the wall should be evalua 
ted by stress linearization. 


Inelastic Solution for One or Two Cycles 


If better accuracy of the results is required, the 
resultants may be obtained by inelastic analysis of one 
or two cycles. The thermal analysis remains the seune. 
Instructions for performing the inelastic stress and strain 
analyses using the MARC program are given in Appendix B. 

The program optimizes the load steps to obtain accurate 
results within the shortest computer time. As in elastic 
euialysis, the bending and shear stress resultants are 
determined simply by linearizing the stresses through the 
thickness of the ligament. The range of plastic strain 
in the hoop direction is determined by adding absolute 
values of strains averaged through the thickness of the 


ligeunent within the entire cycle and taking half of the 


sum as ji strain range. 
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3.6 Therm&lly Induced Stresses 

The h . story of thermally induced stresses for the ana- 
lyzed ligcunent is shown in Figure 8. As explained before, 
the ligaments are subjected to severe thermal straining due 
to the temperacure difference between the outside and inside 
of the thrust chamber. As can be seen in Figure 3, the..e 
is also an essential temperature drop through the wall of 
the ligament. This causes bending since the ligament is 
constrained at the ends. The bending stress arises dv t 
the initial cooling phase, when the rapid chilldown Cc. 
temperature gradients through the ligeunent. Tensile me., 
brane stresses in the ligeunent exceed yield. The stresses 
increase until the temperature difference between the 
ligament and closeout wall reaches a maximxun, after which 
unloading occurs. The stresses become compressive in the 
ligeunent in the residual state when a uniform temperature 
is attained. 

Following ignition, the inner surface of the cylinder 
is subjected to intensive heat flux through the ligament. 

The temperature drop through the ligament results in ther- 
mally induced bending stresses. Compressive membrane 
stresses in the ligament exceed yield. Plastic flow in the 
hoop direction erases the effect of bending and the stresses 
remain essentially uniform through the ligcunent. Unloading 
occurs after the temperature difference between the ligament 
and closeout wall reaches a maximum and the stresses are 
relieved. The bending stresses generated in the ligament 
during unloading are again erased by straining in the hoop 
direction as the entire nozzle cools to a iinifon< temperature. 
The repetitive distribution of stress in subsequent cycles are 
illustrated in Figures 8a and b. 
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a) _ Beginning of Heat Transient Peak of Heat Ttansient 



FIGURE 8 REPETITIVE CYCLE OP THERMALLY 

INDUCED BENDING STRESS IN LIGAMENT 
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The thexnially induced heading though acting only for a 
short portion of the cycle, may enhemce the ratcheting strain. 
Its effect may be assessed by coaqputing the elastic energy of 
the thermally induced bending stresses and correcting the 
hoop strain accordingly. 

The maximum thermal bending stress is given by equation 
(25) , where dT is the temperature drop across the ligament. 
Thus, referring to Figure 9, the elastic energy is: 


E 


el 



(27a) 


_ EH(g A T)^ 
6(1 - v )2 


(27b) 


It is conservatively assumed that all of the avail* 
able elastic energy goes into plastic straining of the 
ligament. 

If Ae” denotes the corresponding plastic hoop strain. 

Pi 

then the plastic «rork is: 



2 SyiCp. H 


Equating (27a) and (28) : 


i28) 


Ac " = ^(** ^ *^*) ^ (29) 

Pi 12(l-v)2Sy 


where 
and S 


is the average of the absolute values of S 

ft 

^min 


y 


max 


This must be added to the plastic strain range due to 

differential thermal expansion, to obtain the plastic 

strain ranac in the hoop direc m, • 

^1 
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FIGURE 9 THERMAL BENDING IN LIGAMENT 
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4.0 DEVELOPMENT OF FAILURE MODE EVALUATION CRITERION 


The ligaments are subjected to incremental plastic deforma- 
tions during each firing cycle of the thrust chamber. The 
geometry of these ligaments changes as the incremental strains 
accumulate. They are subjected to incremental inward bending 
emd simultaneously to a progressive thinning near the center 
of the ligament. Both fatigue damage and the tensile stability 
of the material are affected by local geometry changes. 

Plastic tensile instability occurs when the incremental 
strain hardening of the deforming metal is less than the incre- 
mental increase of true stress due to local thinning. Once 
conditions of tensile instability are reached, further stretch- 
ing occurs as minimum ligament section nec’:ing. The pressure in 
the chamber causes ultimate failure of the necked ligament. 

4.1 Plastic Insted)ility - Necking 

The plastic tensile instability of the ligament in the 
displacement controlled thermal cyclic strain field is 
cmalyzed first by considering the ligament as a biaxially 
loaded shell subjected to monotonic tensile straining in 
the hoop and axial directions. 

For a biaxially stretched sheet, the Mises yield con- 
dition is: 

(30a) 

where ? is the effective stress and and O 2 ate principal 
stresses in the hoop and axial directions, respectively. 

®1 

Denoting the stress ratxo — by a, ernation (30a) can 
be written as: 

a = Oj^d - o + a^)^ (30b) 
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In the present analysis elastic strains will be neglect- 
ed and the Levy-Mises flow rule equations assximed to apply. 
The relation bet%»een the strain increments can then be 
%nritten as [6] : 


de 2 de 2 


dE3 de 


20 ^- 02 20^-02 


Ola) 


tidiere 


dc 




^ +de 2 ^ +de^^j 


(31b) 


In the above, d7 is the increment of effective strain and 
dEj, de 2 aiid de^ are the strain increments in the hoop, 
axial and radial (thickness) directions, respectively. 
Equation Ola) can be written in terms of a as: 


de 


de. 


-de 


de 


2-a 2a - 1 l+o 2(l-a + a2)% 


Ole) 


For proportional loading the stress ratio remains con- 
stant and equation Ole) can be integrated to give the 
total strains: 


- e . 


2-0 2o-l 1 + 0 2(1 -o+o2)^ 


Old) 


Instability occurs when the hoop force F reaches 
a inaxirotam, i.e.: 


<3F-d(2Hao2) *0 


02 ) 


where 2H denotes the thickness of the ligament and a its 
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axial length. Equation (32) is differentiated to give: 



da dH 
a“ ” 'h* 

~ - de_ ~ de 
2 3 


(33) 


Frcxn the incompressibility condition of plastically 
deromied material: 


or 


de + de + de = 0 
12 3 




(34a) 

(34b) 


For proportional loading: 
do j do 



Thus, equation (33) becomes: 



(35) 


(36) 


Using the flow law equation (31c), the condition for 
instability is given by: 


^ - (2 - «) 0 
dl 2(3 - a + 


If the stress-strain law is given by: 

_ _ n 

o = Ae 


(38) 
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t^ere A and n are material constants » then on differenti- 
ating : 

do _ n o 
dc e 

do _ q 

dl~ ^ (39) 


%rtiere z - e/n is the subtangent in Figure 10. 

From equations (37) and (39) : 

— 2(1 — a + ^ 

^ ' (2 -a) 


(40) 


and the critical effective strain is: 


- ^ 2nd -g 


cr 


(2 -g) 


(41) 


The critical strain in the minimum ligament section from 
equation (31d) is then: 

c, = n (42) 

cr 

For cyclic loading, the compressive straining of a 
ligament at the beginning of the loading cycle is followed 
by tensile strains. In each loading cycle the ligament 
accumulates plastic strain. The net increments in the 
hoop and radial directions occur at the end of the cycle 
when the axial strain is zero. Thus, the limits of 
material stability C 2 ui be approximated by the analysis 
for plane strain conditions. The plastic strains in the 
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hoop direction and in the radial direction de^ satisfy 

the condition of incompressibility. Since de^ = Os 


de 


1 



Now, 



(43) 


(44) 


^ere 2H denotes the original ligament thickness and t 
the ligament variable thickness. Then using equations (42), 
(43) and (44) , the critical thickness at instability is given 
by: 

t = 2He cr 


t 


cr 


2H e"“ 


(46) 


4.2 Strain Range - Fatigue 

For fatigue calculations, the maximum local hoop 
strain reuige at the minimum ligament section can be obtain- 
ed from the average hoop strain range by integration if the 
geometry of the distorted ligament is known. As mentioned 
earlier, tests indicate that the deformed shape of the 
ligament can be approximated by a linear function as shown 
in Figure 6. Thus, the thickness variation can be written 
as: 

t(x) = ft - t . + t (47) 

V roax roin / i ^ ^min 

where t emd t . are as shown in the figure. Now, the 
max mxn ^ 

average strain in the hoop direction is given by: 
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■Ij: 


(x) dx 


From equilibrium of forces in the hoop direction: 


c ix) t(X, = 0 t ^ 

i * ml in 


Using the yield condition equation (30b) , since a is con- 
stant equation (49) becomes: 

o(x) t(x) = Vin^min 

Substituting the stress-strain law, equation (38): 

F(x)" t(x) = F " t . (51) 

'' ' * min min 

Now ma)cing use of the flow law, equation (31d) for constant a: 

p : (x)” t (x) = r n t . ^52) 

min 

Substituting for e^(x) into equation (48), the relation 
between the average strain and local strain is given by: 


^min I / ^min v 

— JltOTj 


(53a) 


Substituting for t(x) from equation (47) and integrating: 


t r 2 -=^ 

, - e, jgj-n /jnax\ " 

avg min ” ( ^max “ ^min ) \ *^min / 


(53b) 
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or 


1 


= E . 


xnin 


avg 


\ mxn 


n - 1 -I 


/ 


max 


min 


n 


- 1 


(53c) 


As the deformation increases with each cycle, the thick- 
nesses t and t . also change per cycle. This, in turn, 
max mm 

changes et , the local strain in the minimum ligament 
■^min , . . ^ 

section wnicn intensifies with each cycle. 

The axial strain in the minimum ligament section is 
given by: 

e, = a(T. -T ) (54a) 

2min ^ ° 


Also, from the condition of incompressibility: 


^3, 


min 


= “(e 


Imin ^ *’2min^ 


(54b) 


Now the effective strain range in the minimum ligament 
section for entering the fatigue curve is: 


P 


'min 


" f y^^lmin “ ^2niin^^ ^^Imin ” ^3min^^ ^^2min " ^3min 

(55a) 


Substituting equation [54b) into (55a) ; 

Emin® 3 ^^^Imin ” ^2niin^ ^ ^^^Imin ^2njin^ ^^^2ndn^ ^^min^ 



^lmin^2nain 


®2min 


(55b) 


34 



For the linear variation assumed, the thickness t . 

th min 

and t^^^ after the N cycle are given by: 


2HU -»■ w) - N 6 w 
min (a + w) 


(56a) 


2H(t -<• w) ^ t N 6 & w 
max ~ U + w)2 


(56b) 


where 6 is the deformation per cycle. 

Equations (55) and (56) along with the fatigue curve 

can be used to determine the fatigue life. However, since 

ex . changes with each cycle, a numerical procedure is 
min 

necessary for performing the calculations to determine 
fatigue life. A FORTRAN program has been written for 
determining the cycles to failure and is included in 
Appendix C. 


5.0 THRUST CHAMBER LIFE PREDICTIONS 

The results obtained using the analyses described can be 
used to determine the life of the thrust chamber. For small 
ligament distortion it can be conservatively assumed that the 
incremental thinning of the ligament remains constant during 
subsequent cycles. The number of cycles to failure can then 
be bounded from below by considering the number of cycles 
needed to reduce the thickness of the ligament below the 
critical value resulting in tensile instability or by con- 
sidering the fatigue damage of the gradually thinned ligament, 
whichever results in smaller number of cycles. For OFHC 
copper this bounding technique provides a realistic evalua- 
tion of the cycles to failure observed in experiments where 
failure was due to tensile Instability. 
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For NARloy Z the thinning gradually diminishes in con- 
secutive cycles and the fatigue mode of failure prevails. 

This bounding technique provides a conservative evaluation. 
However, the efficiency of the bound is significantly lower 
than for copper. 

The tendency of diminishing increments of accumulated 
strain in consecutive cycles is the result of kinematic hard- 
ening. It has been demonstrated [7], [8], that for such 
material the cyclically loaded structure always achieves the 
condition of alternating plastic straining if there is no limit 
on the hardening capacity and the changes of geometry can be 
ignored. The net increment of plastic ratchet strain vanishes 
after sufficient hardening. The cunount of plastic strain 
accumulation needed to achieve such plastic shakedown can be 
uniquely determined. For real materials with limited harden- 
ing capacity and less vhan ideal kinematic behavior, the 
stable cyclic state occurs at much larger accumulated strains 
than for an ideal kinematically hardening material or may 
never be achieved. Materials more closely obeying the theo- 
retical concept of kinematic hardening can be expected to fail 
due to fatigue since slowed down thinning allows the material 
o be exposed to more fatigue cycles, as in the case of 
NARloy Z. 

This trend is also in general agreement with the results 
of (91. The number of cycles at which thinning stops is 
related to the strain hardening parameter n. Assuming this to 
be the dominating parameter on whic’ the number of cycles to 
thinning depend and hypothesizing a power relationship, the 
following empirical criterion is obtained: 

= 750n^*^^ (57) 

where denotes the nximber of cycles at which thinning stops. 


36 



Note that the results of [9] have been partly used in determin- 
ing the constants for the power relationship since these were 
the c'ly quantifiable numbers available. 

If the number of cycles to reach the critical thickness 
obtained from plastic instability analysis is less than N^ for 
that material, then tensile stability is the limiting failure 
mode. The thrust cheunber life in this case is given by the 
number of cycles to instability. If on the other hand, N^ is 
greater, fatigue is the failure mode and thrust chamber life 
is > iven by N^ plus the cycles obtained from the fatigue curve 
corresponding to the effective strain range which remains con- 
stant from that point on. 

The numerical calculations for the case of OPHC copper and 
NARloy Z using the procedure developed herein are presented in 
Appendix A. 

6.0 FINITE ELEMENT ANALYSIS 

The MARC program was used for the thermal and mechanical 
part of the analysis. The MARC program was adapted to analyze 
the fast thermal and pressure thrust nozzle transient. Opti- 
mized selection of the time ar.u load steps was used for the 
theinnal and elastic-plastic analyses. 

The finite element inelastic analysis confirmed the fact 
that the simplified model includes the essence of the physical 
behavior of the ligament as described herein during the course 
of the development procedure. Moreover, the results of the 
finite element analysis were employed in deriving and cali- 
brating the thermal model of the ligament. A description of 
the model and details of the finite element analysis for the 
case of OFHC copper are given in Appendix B. 
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7.0 COMPARISON OF RESULTS 

The numerical results obtained for OFHC copper, using the 
simplified analysis procedure developed herein, are compared 
with the finite element results in Table 1. The results are 
seen to be in good agreement. The life prediction of 103 cycles 
for OFHC copper and 833 cycles for NARloy Z is also in general 
agreement with experimental results [2] indicating that the 
simplified design procedure can be used for predicting the life 
of the thrust chamber without performing detailed inelastic anal- 
ysis. 

TABLE 1 

Comparison of Analytical vs. Finite Element Results 
Residual Deflection After One Cycle 


Analytical 

Results 

Finite Element 
Results 

Deflection of Pres- 
sure Surface/Cycle 

Deflection of Pres- 
sure Surface/Cycle 

.00029" (.00737 mm) 

.00027" (.00686 mm) 


8.0 CONCLUSIONS 

A simplified design procedure for predicting thrust chamber 
life is developed herein. The method uses a yield surface for 
combined bending and membrane loading to determine the incre- 
mental inward bulging and progressive thinning near the center 
of the ligaments at the inner liner of the thrust chamber. 
Failure analyses indicate that plastic tensile instability is 
the dominant mode of failure for OHFC copper. Both fatigue 
and plastic tensile instability must be analyzed for NARloy Z 
to determine the limiting failure mode. Results of the 
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simplified analyses are shown to compare favorably with those 
obtained from detailed inelastic finite element analyses for 
OFHC copper. 

Further experiments are needed to provide better characteri- 
zation of plastic response of the materials used for thrust 
ch2unbers. 

For longer cycle times, the accxiracy of the evaluation can 
be improved by including the time dependent response of the 
material into the analysis. 
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APPENDIX A 
NUMERICAL EXAHPL s 


Exanqple l 

Consider the ligament for OFHC copper. The geometry is as 
shown in Figure A-1. 

For a unit width of the ligament Ap = 547 Ib/in. (95794.47 N/m) 
The average yield stress Sy = 9000 psi. (62 MPa) 

Then from elastic solution for a clamped beam (eqs. (19) and (20)) 


Moment M 

Shear Force S 

From eq. (3a): 

Yield Moment M^ 

o 

From eq. (2a): m 

IP 

m 



Ei - 

2 


px 


S. 


JL 

M_ 


-547 . (.0664)2 




9000 


fi £ 

|_6 .06 


0664 (.0664) 2 


-0.0729 + 6.5888X - 99.229x2 


(A-1) 


From eq. (10b) : 


s 


s 



m 


547 . .0664 


- 547x 


9000 


= .1153 - 3.473X 


(A-2) 
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From eq. ( 6 b) : 


k = 1 - 

= (.9867 + .8009X - 12.062x2)** (A-3) 


From eq. ( 6 a) the yield siirface: 

fj2 

F =m+2^-k=0 

n2 = r 2 - km (A-4a) 


Substituting eqs. (A-1) and (A-3) : 

n = ^.9867 + .8009X - 12.062x2 

- (.9867 + .8009X - 12.062x2)** 

(-.0729 + 6.5888X - 99.229x2) 1 (A-4b) 

From eq. (21), taking a. = o - a and (Sy . )/2=S , the 

2 . o ->^max ^mxn ^ 

plastic strain range in the hoop direction is: 

Ae* = fa{(T.-T) -(T.-T) . } 

Pj^ L * ^ o'inax X o^mxnf EJ 

(A- 5a) 


From thermal analysis, (T^ ' J =780"F ( 433 “K) 

Then tor a = 9.5x10“^ in./in./®F, ( 17,1 x 10~® ®C“^) and 

E = 17x10® psi (117215 MPa): 


lie 


Pi 



10“® 


780 


2 ♦ 9000 
17 • 10® 


Ae • = 0.C0635 
Pi 
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From eq. (29) : 


Ac • = Ea2(AT)2 

12(1 - v)2Sy 

The temperature drop across the ligament AT ' 200*F: 

^ 17 • 10^ • (9,5) ^(10’^) ^(200)^ 
'Pi 12(1 - 0-3)2 . 9000 


(A-5b) 


Ae " = .00116 
P, 


.*. Ae = Ae ' + Ae • = .00751 

Pi Pi Pi 

From eqs. (6b) amd (15); 



(A-6) 


(A-7) 


For the present problem, m and s are small ccmipared to n. 
Thus it is found from equation (A-4a) that n ' k. Eq. (A-7) 
then sin^lifies to: 


K » 



(A-8) 


K = i = 2 (.00751) 
~ R ~ .0175 


where the radius R is assumed for the length shown in 
Figure A-2. 

Thus R = 1.1651" (29.59 mm) 
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2R(1 > cos d) 


Deflection 6 ^= 

ij 

= 2(1.1651)^1 - { (1.1651) ^ - (.0166)2) /i.iesij 
=. 0002365" (.006 mm) 

From eq. (17), since the contribution of m and s is small, 
the product term can be neglected and n taken = 1 as 

^1 - s2 

seen from Fig. 5b. 

Thus: Y = 4s|ACp j (A-9) 

Substituting from eqs. (A-2) and (A-6) : 

V = 4(0.1153 - 3.473X) (.00751) 

Y = .003464 - 0.104328X (A-10) 

Integrating and determining the deflection at x = t/2: 

= .003464X - 0.104328x2/2 
= C. 0000575" (.0015 mm) 

.*. Total 6 ” ^2 ~ -000294" (.0075 mm) 

From eq. (18) , thinning after N cycles: 

^ _ 0.000294, • 0.05N 

N ” (0.0664 + 0.05) 

tjj = 0. 000126 3N inches (A-11) 

Material Instability 

The strain hardening parameter n in the stress-strain 
law is approximately given by LeRC as: 
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where denotes the ultimate strength and the yield 
strength . 


Towards the end of the cycle, the value of and 
Sy are approximately 46 ksi (317 MPa) and 9 ksi (62 MPa), 
respectively. Thus: 


n 


0.2 



0.6 


= 0.467 


From Eg. (46), critical thickness for n = 0.467 
ter = (0.035)e”°-‘*®'^ 


tej. = 0.02194" (.557 mm) 


Thinning for 

Instability = (.035 - .02194) 


= .01306" (0.332 mm) (A-12) 

Equating (A-11) and (A-12) , the number of cycles to insta- 
bility are: 


N = 103 

From equation (57) 

750 (0.467)^*^^ 

= 289 cycles 

Thus thinning continues for 289 cycles. 

Since N < N^, the failure mode is material instability 
and thrust chamber life equals 103 cycles. 

Fatigue 

To demonstrate that fatigue results in a higher number 
of cycles, the procedure described in Section 4.2 was used to 
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determine the fatigue life of the OFHC ligament. The 
fatigue life was based on the fatigue curve for OFHC copper 
shown in Figure A-3. The calculations were performed using 
the FORTRAN program listed in Appendix C. The fatigue life 
was determined to be 132 cycles. 
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total strain range (PERCENT) 



FIGURE A-3 TYPICAL LOW-CYCLE FATIGUE LIFE OF OFHC 
COPPER ANNEALED CONDITION 
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Example II ; 


Consider now an example of NARloy Z. The geometry is the 
seune as shown in Figure A-1 

Also, Ap = 547 Ib/in. (95794.47 N/m) 

The average yield stress Sy = 30,000 psi (207 MPa) 


As before. 


Pi 


from Eq. (21) 

a j(T. - T ) 

( 1 o max 


(T -T ) . 

1 o min 



(A-13) 


Since the thermal conductivity for NARloy Z is approximately the 
same as that of OFHC copper, temperatures Tj. an T^ are approxi- 
mately the s£une. 

•Hien for a = 9.5 x lO”*' in./in./“F, (17.1x10"® °C"*) and 
E = 18x10® psi (124110 MPa): 

Ae' - 9.5 • 10"® • 780 - - ■ 

Pi 18 • 10® 

Ae' * .004077 

Pi 


From Eq. (29) 


Ae" = 18 • (9«5)^ (10"^)^ (200)^ 

Pi 12(1-0.3)^ • 30000 

Ae" * .000368 

Pi 

Ae„ * Ae' + Ae" = .004445 
Pi Pi Pi 


Substituting in : 

K = 

H 

^ _ 1 _ 2(. 004445) 

R “ .0175 

R = 1.9685" (50 mm) 
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Deflection 6 ^ • 2 R(l-cosa) » .00014" (.0035 mm) 

From Eq. (10b) and (20) 


8 « 


s - 


HS, 




/547 • 0.0664 


30000 


I 


- 547x 


= 0.0346 - 1.042X 
Substituting in: 


Y » 4s (Ae^ ) 

Y = 4(*0346 - 1.042X) (.004445) 

« 0.000615 - 0.018527X 

Integrating and determining the deflection at x*Jt/2 
^2 = 0.00001" (.000254 mm) 

Total 6 “ 6 i + 62 * 0.00015" (.00381 mm) 

Prom Eq. (18), thinning after N cycles: 

♦, ~ 0.00015 • 0.05N 

^ - (0.0664 + 0.05) 

“ 0.0000644 N inches 


(A-14) 


(A-15) 


(A-16) 


Material I itabilitv 


For S- » 55 ksi (379 MPa) and S = 30 ksi (207 MPa) 


u 

n 

n 


0.2 ( 

0.18 


55-30\ 
30 ' 


0-.6 


From Eq. (46) 


t 


cr 


(.035)e‘°*^® 
0.0292" (.742 mm) 


Thinning for instability 


0.035 - 0.0292 
0.0058" (.147 nun) 


(A-17) 
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Equating (A-16) and(A-17), the number of cycles to instability 
N = 90 

From Eq. (57) 

= 750 ( 0 . 18 )^*^^ 

= 38 

Thus thinning stops after 88 cycles. 

Since N > material instability is not the failure mode and 
the thrust chamber life must be determined from fatigue calcu- 
lations . 

Fatigue 

The procedure of Section 4.2 was used for determining the 
fatigue life of NARloy Z. The fatioue life was based on the 
fatigue curve of Figure A-4. The FORTRAN Program listed in 
Appendix C was used for performing the necessary calulations . 

The input and output obtained are attached in this Appendix. 

The notation for the listed results is as follows: 

NCYC - Number of cycles 

TMIN - t_. , minimum ligament thickness - in. 

mxn 

TMAX - t^ , maximum ligament thickness - in. 
max 

EEQVT - effective strain range % 

NF - Nun^er of cycles to failure from fatigue curve 
USEAGE- Useage factor 

Thus if the ligament were to continue thinning, the fatigue life 
would be 272 cycles. However since thinning stops after 88 
cycles, the remaining fatigue life is based on the effective 
strain range of 2.2416% which now remains constant and equals 
747 cycles. 

Thus cycles to failure * 88 + 747 

= 835 
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TOTAL STRAIN RANGE (PERCENT! 



53 


constants input 

N« ,16000 
delta* .15000E-OS 
E1AV6* ,74100X 

H« ,17500E-Ot 
L* ,66400E«01 
Ws ,50000E-01 


NARloy-Z 

n (material strain hardening exponents) 

S (Deflection per cycle, in.) 

^lavg (Average hoop strain in ligament, aAT) = c2 

min. 

half ligament thickness, in. 
t (width of ligament in hoop direction, in.) 
w (width of rib, in.) 


FATIGUE CURVE 
CYCLES STRAIN RANGE 

10.000 16,000 

100.00 5,4000 

4000.0 ^,0000 

.'S0000E«06 ,10000 
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at!«OLlS fUt> r*CH CVCI C |NCTC« TmIN« rN«V,M QvT,Nf #USMUf 


♦ 

. I4436E-U1 

.'ssoiTr.oi 

1.4090 


.57101t*0l 

i 

. 34aTl«-0t 

.iS074f.Ol 

1,4940 


,1 14446*02 

3 

, 34a07t«ut 

.isnoF.oi 

1,50U1 

iJii.v 

•172766*02 

41 

, \4T4?(*U1 

,5Sl47f*Ol 

1.5092 

IJ07.S 

,231326*02 

s 

. 3467ae*0i 

,351446 .01 

1*5124 

I*. 1.0 

•290396*02 

6 

,34M31-u1 

.3S22IF-01 

1,5100 

I.TB.o 

,349966*02 

7 

. 34S491^^'» 

,i5257F*0| 

1,5249 

(■fell . s 

,4lOU4t*V2 

A 

. 344ase«oi 

, JS294F-0I 

1.5312 

(•50. « 

, 4 TObSt *02 

9 

, 34424>i»0l 

. 353S1F-01 

1.5376 

IfeSS.B 

,531746*02 

10 

. |4 9S6(-01 

.3S360r*0l 

1.5441 

IfeJl.fe 

,593456*02 

1 t 

.3429U*U1 

.3S404F-01 

1,5599 

1*07. i 

,655966*02 

12 

. 3422T£*01 

;39«4|F*0| 

1,5571 

14,1. fe 

•710416*02 

1 3 

. 341626*01 

,iS476F*0| 

1,5037 

Ii7fe.i 

•701726*02 

141 

^ 140941 *01 

.3S515F-01 

1,5704 
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APPENDIX B 


FINITE ELEMENT ANALYSIS 


INTRODUCTION 

The MARC general purpose finite element computer program 
was adopted to study the response of the experimental thrust 
chamber, during its firing tests. Each test consisted of a 
chilldown cycle and a niaraber of repetitive firing cycles. 

Finite element heat transfer and structural analyses were 
performed on the model. In the heat transfer analysis, the 
temperature distributions vs. time were obtained for a typical 
chilldown and firing cycle. In the structural analysis, the 
resultant temperature disbributions and the cyclic pressure 
loads were used to determine the linear and non-linear response 
of the structure. Optimized selection of the time and load 
steps was used for the thermal and elastic-plastic analyses. 


THERMAL ANALYSIS 

1) Model Description 

The finite element mesh used in the analysis is 
shown in Figure B-1. This model was established based 
on the geometric symmetry planes and loading conditions. 
It consists of thirty-five elements and 138 nodal points. 
All elements are the eight-node planar bi-quadratic 
quadrilateral type (element 41 in MARC element library) . 
Temperature and position (coordinates) within these 
elements are interpolated from eight sets of nodal 
values, the four corner nodes and four mid-side nodes. 

The interpolation function is such that each edge has a 
parabolic variation along itself. 
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2 ) Material Pro‘perties 

The properties used to characterize the response of 
the nozzle in the thermal analysis were obtained from 
Reference 3. They were thermal conductivity and 
specific heat as a f\inction of temperature of the thrust 
chamber material as shown in Figures B-2 and B-3, 
respectively. 

3) Load Cycles 

The thermal analysis consisted of a chilldown and a 
firing cycle which lasted 1.5 and 3.5 seconds, respec- 
tively. At the beginning of the chilldown period, the 
room temperature value of 294“K (530®R) was assumed to be the 
uniform initial temperature of the entire thrust chamber. 

The thermal operating conditions (heat transfer coefficient) 
were defined in terms of thermal boundary conditions on the 
cooling channel surfaces (cold side) and the combustion gas 
surface (hot side) . These boundary conditions for the chill- 
down and a firing cycle (totally five seconds) are plotted vs. 
time in Figures B-4a through B-4D. To improve the convergence 
process, the temperature of the coolant flow at the beginning 
of the chilldown period was assumed to be the same as the 
initial model temperature 294“K (530®R) , and then it was 
suddenly dropped to 29 ®K (53 ®R) in 0.1 second. For more 
accuracy in results of the heat transfer analysis, the time 
steps were controlled by the program to limit the maximum 
temperature change at any node to 28*K (50*R) . 

4) Results of Heat Transfer Analysis 

The heat transfer analysis of the model was completed 
in 107 iteration cycles, 25 iteration cycles for the chill- 
down and 82 iteration cycles for a firing cycle. Results 
show that the temperature of the entire model drops from 
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SPECIFIC HEAT J/ cm V“K (BTU/IN V°R)x10-^ CONDUCT! VI TV W/cm/’R (BTU/IN/S/"'R)xlO"’ 


ORIGINAL ® 

OF POOR QUALITY 




FIGURE B-3 MATERIAL SPECIFIC HEAT VS. TEMPERATURE 
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FIGURE B-4 THE TIME HISTORIES OF HEAT TRANSFER COEFFICIENTS 
AND TEMPERATURES FOR A CHILLDOWN AND FIRING CYCLE 
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294®K (SaO^R) at t = 0 to less that 30“K (54«R) at 
t = 0.65 second and to a uniform temperature of 29®K 
(53®R) at t = 1.13 seconds. From t = 1.13 to t = 1.5 
seconds the nozzle remains at 294®K (530®R) temperat re. 
From t = 1 . 5 seconds which is the start of the firing 
cycle, the temperature of the model begins to rise until 
the temperature of the hot surface reaches 805®K (1450®R) 
at t = 3.2 seconds. At this point the temperature of 
the model begins to drop until the entire model reaches 
the temperature of 29“K (53“R) at t = 5.5 seconds. 

To get a picture of the temperature distribution, 
plots of continuous isotherms at t = 3.2 seconds are 
reproduced in Figure B-5. The temperature distribution 
of the centerline of the model is also plotted at 
different time intervals in a firing cycle in Figure B-6. 
It appears from the results that the maximum temperature 
difference between the hot and cold surfaces is 114 ®K 
(206 ®R) and the maximum temperature difference between 
the hot surface and the outer surface of the model is 
^25®K (946“R) during the time interval of t = 1.9 to 
t = 3.2 sec. The temperature histories of a few selected 
points are also plotted for a chilldown and firing cycle 
in Figure B-7. Since the temperature of the entire thrust 
chamber approached 29 ®K (53®R) at the end of the first 
firing cycle, it was concluded that the thermal analysis 
of this cycle without the chilldown could be used to 
simulate the remaining cycles. 
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TEMP. DIST. RUN 42 


FIGURE B-5 ISOTHERM PLOTS AT t=3.2 SECONDS 
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FIGURE B-6 MIDPLANE TEMPERATURE HISTORY 
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STRUCTURAL ANALYSIS 


1) Model Description 

The finite element mesh used in the structural 
analysis is the same as the one used for heat transfer 
analysis. Eight-node isoparametric distorted quadri- 
lateral generalized plane strain elements are used for 
the structural analysis. Displacement and position 
(coordinates) within the element are interpolated from 
eight sets of nodal values. The four corners and four 
mid-side nodes. The interpolation function is such that 
each edge has parabolic variation along itself. 

2) Material Properties 

The properties used to characterize the response of 
the thrust chamber in the structural analysis were taken 
from Ref. 3. They were time-dependent thermal expansion. 
Modulus of Elasticity, and stress-strain curves (work 
hardening) . 

3) Loading Cycles 

The temperature distribution from the heat transfer 
analysis and the corresponding cyclic pressure load were 
supplied to the structural analysis program. The cyclic 
pressure load as a function of time is plotted in Figure 
B-8. The structural analysis was performed for a chill- 
down and two consecutive firing cycles. Due to the non- 
linear response of the thrust chcunber material, iteration 
cycles were controlled to limit the maximum temperature 
change input to the structural analysis to 13.9*C (25"Fjat 
any node, loth isotropic and kinematic work hardening 
were used xn the analysis. 
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FIGURE B-8 CYCLIC PRESSURE LOAD VS. TIMc 
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4) Results of Structural Analysis 

The structural analysis of the model was completed 
in tliirty-four (34) iteration cycles simulating the 
chilldown and 130 iteration cycles simulating each firing 
cycle. Some plots of typical stress and strain distri- 
butions are reproduced in Figures B-9 through B-14 . 
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CURVE NO. (PSi) 



FIGURE B-9 STRESS INTENSITY (VON MISES) DISTRIBUTION FOR THE 
MAXIMUM PRESSURE DIFFERENCE (ELASTIC ANALYSIS) 
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CURVE NO. 


FIGURE B- 


(psi) 






TIME=1 .7 


STRESS INTENSITY (VON MISES) DISTRIBUTION AT TIME =0.2 SECONDS AFTER 
THE BEGINNING OF THE FIRST FIRING CYCLE (ELASTIC ANALYSIS) 
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FIGURE B- 
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TIME=1 -7 


STRAIN DISTRIBUTION IN THE RADIAL DIRECTION AT TIME = 0.2 SECONDS 
AFTER THE BEGINNING OF THE FIRST FIRING CYCLE (ELASTIC ANALYSIS) 
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ORVE NO. (psi) 



FIGURE B-12 STRESS INTENSITY (YON NISES) DISTRIBUTION AT 
THE END OF THE SECOND FIRING CVaE 


original is 
OF POOR QUAUTY. 
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OnOMALPAOEB 

OF POOR QUALITY 

CURVE NO. 




I i TIME---8.5 

FIGURE B-13 TOTAL EQUIVALENT PLASTIC STRAIN DISTRIBUTION AT 
THE END OF THE SECOND FIRING CVaE 
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CURVE NO. 



FIGURE B-14 PLASTIC STRAIN DISTRIBUTION IN HOOP DIRECTION 
AT THE END OF THE SECOND FIRING CYCLE 
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APPENDIX C 


FORTRAN Program 
> for 

Fatigiie Calculation 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


PliOCRAN CYClM 74/7R OPT«| 
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the format for READING THERE CONRTANtb IR bf|o 
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IR no 

CTCIEacvclES for Each data PT CtFloi rpt card if NPTR GT r 
STRAINaRtRAIn range for each PT (RFtoi RPT CARO IF NPTrIgtIr 

READ too, NPTR 
READ 200, (CYC(.E(I},2 p|,NPT 8) 

READ 200, (RTRAtNmflAl,NPTR) 

PRINT SOO/ N,OClTA,EIAVG#M,l,N 

DO S lAl^NPta 

S PRINT AOO* CVClC(n,RTRA|NCt) 

calculate CTCLER to failure 

PRINT 700 

useageao.o 

NCTCaO 
Ca(N-1)/n 
10 NCYCANCVC^t 

TNXNA(2*M*(LAU)«NCVCA0ELTAftH|/(L0M) 
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NE9ULTS 

POINT «oo, Ncvc-i.usiose-i./NP.Ncvc.unoce 

STOP 


program cycles 74/70 OPTa| 


FTN «,4#RS2 


40 


45 
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lAu FORHATCIIO) 

200 FORMATdFSOtO^ 

300 FORNATClHO#ATHC STRAIN RANGE OF AfGtSjS#* FOR THE *.16#* CVCLf lA 

I not on the fatidue curve SUPPLIEOa) ' 

350 FORHATClMOfATNE USIAGC FACTOR IS AiFaJoJ 
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APPENDIX D 


SYMBOLS 

a 

A 

D 

E 



F 

h 

2H 

k 

K 

£ 

m 

M 

M 

o 

n 

N 

P 

s 

S 
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axial length 

surface area, also equation arbitrary positive 
scalar 

rate of dissipation 
modulus of elasticity 
elastic energy 
plastic work 
yield surface 

convective heat transfer coefficient 
thickness of ligament 

- s^ ; also thermal conductivity 
curvature 

width of ligament in hoop direction 
generalized bending stress variable 
bending moment in ligament 
yield bending moment 

generalized hoop stress varicible; also strain 
hardening exponent in stress strain law 

hoop force in ligeunent; also number of cycles 

yield hoop force 

pressure 

generalized shear stress variable 
shear force in ligament 



SYMBOLS ” continued 


So 

— 

yield shear stress 


- 

average yield stre- ■ in tension 

s„ 

- 

ligament yield strength for minim’im o(T. - T ) 

■^min 

S„ 

- 

ligament yield strength for maximum a(T. -T ) 

-^max 

t 

— 

ligament thickness 


- 

thinning after N cycles 

t 

max 

- 

Fiaximtua ligament thickness 

t . 
min 

- 

minimum liqcT'ent thickness 

T 

- 

temperat’ 


- 

average temperature of ligeiment 

T 

o 

- 

average temperature of closeout «rall 

V 

- 

width ■Jf rib 

z 

- 

F/n, subtangent in Figure 10 


a 


coefficient of thermal expansion; also stress 
ratio 

Y 

- 

shear strain 

6 

- 

deflection per cycle 

h 

- 

deflection due to moment 


- 

deflection due to shear 

Ap 

- 

pressure difference between coolant pressure 
combustion gas pressure 

U) 

<1 

- 

total strain range in hoop direction 

Ae 

Pi 

- 

plastx - strain range in hoop direction 
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SYMBOLS - continued 


Ae' 

Pi 

Ae" 

Pi 

e 

e. 


avg 

'^min 


mxn 


min 


£ 

F 


cr 


^min 


G 

X 

V 


a 

4 » 


plastic strain range in hoop direction due to 
differential thermal expansion 

correction to plastic strain range in hoop 
direction due to thermally induced bending 

hoop strain 

hoop otrain 

axial strain 

radial strain 

average hoop strain in ligaunent 
hoop strain in minimum ligcunent section 
axial strain in minimiun ligaunent section 
radial strain in minimum ligeiment section 
effective strain 
critical effective strain 

effective strain in minimum ligament section 

generali'^ed bending strain variable 

generalized hoop strain variable 

Poisson's ratio 

hoop stress 

axial stress 

bending stress 

effective stress 

generalized shear strain variable 
dot above symbol denotes rate 
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